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ABSTRACT 

Scattered light echoes from variable and cataclysmic stars offer one of the most effective 
means to probe the structure and composition of circumtcUar and interstellar media. I build a 
simple model of light-echo surface brightness by considering the source spectrum, and the dust 
density, geometry, and scattering efficiency. I use this model to investigate whether echoes should 
be observable around short and long-period giants, cataclysmic variables, and supernovae. Only 
supernovae are expected to illuminate material on both circumstellar and interstellar scales. Giant 
and post-AGB stars (e.g. Cepheids and Miras) with high mass-loss rates can produce observable 
echoes within their circumstellar envelopes. Echoes from novae and dwarf novae are probably 
detectable only from nearby material, and only in unusually dense gas. I present characteristic 
exposure times to image such echoes over a wide range of wavelengths for ground-based and Hubble 
Space Telescope observations. I apply these results to analyze the dust properties of the recently- 
reported echoes around SN 1993J, finding the dust in M81 to have a grain-size distribution 
and chemical composition consistent with Galactic dust. Optimal observing strategies for echo 
detection are also discussed. 

Subject headings: circumstellar matter — ISM: dust — stars: supernova (SN 1993J) — stars: 
variables - techniques: photometric 



INTRODUCTION 



Since Couderc (1939) proposed the theory of "aureoles lumineuses" to explain the otherwise superluminal 
sources surrounding Nova Persei 1901 (Ritchey 1901), scattered light echoes have offered one of the most 
effective means to probe circumstellar and interstellar structure. When the light pulse from some varying 
object is scattered by dust into the line of sight, an observable echo is produced, provided the pulse is 
sufficiently luminous and the dust sufficiently dense. An echo observed a given time after the pulse must lie 
on the locus of points equidistant in light travel from the source and observer, i.e. an ellipsoid with known 
foci. This simple geometry directly yields the three-dimensional position of an echo, uncertain only by the 
assumed distance to the source. 

A great challenge to any observer is determining the distances of objects of interest along the line of 
sight. Whereas this normally requires a variety of positional and absorption- velocity information as well as a 
geometric model, the relative line-of-sight positions of light echoes are unambiguous, and can be determined 
with as little as a single image. Given a reliable distance determination to the source, the absolute position 
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of an echo and the precise three-dimensional structure of the surrounding dust are known. Conversely, light 
echoes can be employed to determine the distance to structures with independently-known geometry. Finally, 
as the spectrum of an echo is determined by the scattering properties of the dust, light echoes constrain the 
grain-size distribution, density, and composition of dust in the circumstellar (CSM) and interstellar medium 
(ISM), provided the input spectrum is known. 

Light echoes have been considered in connection with gamma-ray bursts (Reichart 2001; Band & Hart- 
mann 1992) and active galactic nuclei (de Diego et al. 1998), as distance indicators to galaxy clusters (Milgrom 
1987; Katz 1987; Braun & Milgrom 1989), and for probing the baryon density of the universe with high-z 
sources (Sholomitskii 1997). If nearly-90° scattering occurs, the positions of maximally-polarized echoes 
yield geometric distance measurements to the echo source (Sparks 1994, 1996). The influence of echoes on 
the light curves of type II supernovae (SNe) has been discussed by e.g., Di Carlo et al. (2002), Roscherr & 
Schaefer (2000), Filippenko et al. (1995), Chugai (1992), and Mackey (1987). The detection of echoes has 
been investigated around SNe (Wright 1980; Chevalier 1986; Schaefer 1987a,b; Emmering & Chevalier 1989; 
Sparks 1994; Maslov 2000), novae (Gehrz 1988), and flaring stars (Argyle 1974; Bromley 1992; Gaidos 1994). 

In spite of the wealth of information contained in echoes, light echo astronomy is still a nascent fleld. 
This is due in part to the relative difficulty of directly detecting them. Foremost, the brightest echo-producing 
objects (SNe) are also the rarest, and only a limited number of recent candidates are close enough to produce 
resolved echoes. Targeted searches around supernovae have been conducted by van den Bergh (1965a,b, 1966) 
and Boffi. Sparks. & Macchetto (1999). The latter search found 16 echo candidates out of a study of 65 
historic SNe in 38 galaxies. Although no conclusive detections were made, echoes could be discovered around 
these candidates in follow-up observations. Novae are the next likely class of targets, and of the searches 
conducted by van den Bergh (1977) and Schaefer (1988), no echoes were detected. To date, echoes have 
actually been resolved around only a handful of objects. These include Nova (N) Per 1901 (Couderc 1939; 
Swenson 1948), N 1975 Cyg (Bode & Evans 1985), N V838 Mon (Munari et al. 2002; Bond et al. 2003), 
SN 1987A (see below), SN 1991T (Sparks 1994; Sparks et al. 1999), SN 1993 J (Liu et al. 2002; Sugerman & 
Crotts 2002), and SN 1998bu (Cappellaro et al. 2001). Variability in the circumstellar material around the 
Cepheid variable RS Pup (Havlen 1972; Mayes, Evans & Bode 1985) and the biplolar nebula OH 231.8-1-4.2 
(the Rotten Egg nebula; Kastner et al. 1992) may also be interpreted as light echoes from the central variable 
star. 

SN 1987A in the Large Magellanic Cloud (LMC) is the nearest known SN in the last 400 years. Due to its 
proximity, echoes around SN 1987A are the most widely discussed and extensively studied. It is beyond the 
scope of this article to give a full accounting of the references-there are at least 15 publications discussing 
the detection of echoes, not including regularly cited circulars and conference proceedings. These echoes 
have been used to map the surrounding circumstellar (Crotts & Kunkel 1991; Crotts, Kunkel & Heathcote 
1995; Sugerman et al., in preparation) and interstellar (Xu, Crotts & Kunkel 1994, 1995; Xu & Crotts 1999) 
material, thereby probing the progenitor's mass-loss history, its location within the LMC, as well as the 
structure and history of the associated stars and gas. 

The low detection rate from intentional searches and small number of known echoes suggests that they 
are especially diflicult to image, and may have discouraged interest in recent years. One may also speculate 
whether limitations in resolution and detector sensitivity, as well as unoptimized observing strategies, have 
contributed to these disappointing results and to the limited application of this technique. Under what 
circumstances can one hope to succesfuUy image light echoes? Are there better observational techniques or 
more likely candidates to study? How can one maximize the information content of those observed? 
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In this article, I discuss the observabihty of scattered-hght echoes around a variety of variable sources, 
and present an optimized observing method, which has already been proven successful. Formalisms for dust 
scattering are developed in §2, in which I build the components for predicting light-echo surface brightness 
in §3. Dust scattering properties are presented in §3.1, candidate variable stars for producing echoes in 
§3.2, and input constraints in §3.3. These components arc used in §3.4 to compute the exposure times to 
image light echoes using ground-based telescopes, and the Advanced Camera for Surveys (ACS) and Space 
Telescope Imaging Spectrograph (STIS) aboard Hubble Space Telescope (HST). I apply the model to echoes 
recently discovered around SN 1993J in §4, and discuss an optimal observing strategy in the conclusions 
(§5). 



2. FORMALISM 

2.1. Light Echo Surface Brightness 

Consider the geometry shown in Figure 1 (adapted from Xu ct al. 1994), in which an object O at 
distance D from the observer emits a pulse of light of duration r at time i. This light echoes off dust located 
at position r from O, and reaches the observer at time t. This formalism was first developed to discuss 
light echoes from supernovae (Emmering & Chevalier 1988, 1989), for which i is measured from the time of 
explosion, and t is measured from the arrival of first visual light. For the more general case of any variable 
source, one can consider i to be the time at which maximum light was emitted from the source, and t is the 
time an echo is observed as measured from the arrival of maximum light. 

Let us treat the dust as a thin planar sheet of thickness Az; this approximation is only for computational 
convenience but will not effect the generality of the result. The light scattered at r lies on an ellipsoid with 
O and the observer at the foci. Since D is generally much greater than r, the echo depth z along the 
line-of-sight can be approximated by the parabola (Couderc 1939) 

where p = rsin^ is the distance of the echo from the line-of-sight in the plane of the sky, and 9 is the 
scattering angle. Since the light pulse and dust have finite widths, the echo will have a width dp (§2.2). 

The fiux scattered off one dust grain of radius a at position r is (Chevalier 1986; Emmering & Chevalier 
1989) 

rw /A ^ -^ - Qsc(A,a)(TgX(A,f)^(M,A,a) 

where Qsc is the grain scattering efficiency; ag = Tia^ is the grain cross section; L{X,i) is the luminosity at 
A and i; and $ is the scattering phase function for a given p = cos0 = 1 — c{t — i)/r. I adopt the Henyey & 
Greenstein (1941) phase function 

$(M, A, a) = l:i9V^ 

[l + 9\X,a)-2giX,a)pf' 

with g{X, a) measuring the degree of forward scattering for a given grain. The total flux Fsc from a single 
scattering is found by multiplying equation (2) by the dust density nrf(r, a) and integrating over the scattering 
volume and all grain sizes. I restrict this work to single-scattering events, and refer the reader to Chevalier 
(1986) for a discussion of multiple-scatterings, which are only of significance when the optical depth of the 
scattering material > 0.3. 
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Proceeding as in Xu et al. (1994), for a thin sheet, changes in z can be ignored over the scattering 
volume. First, note that the flux F{X,i) = L{X,i)/^-iTD'^, then 

Fsc{X,t) = j Qsc{\a)agF{\t)^{p,\a)'^^^^^dad^r. (4) 

Solving n for t and substituting r = \/ (P- -\- z^. 



and for a given time t 



c c 

dt=-'-^. (5) 
cr 

Consider the spatial part of equation (4) by transforming the integral to cylindrical polar coordinates 

nd[r,a) ^^^^ pdp dzd(p. 
Substituting equation (5) and letting dz = Az and F{X) = J F{X, t)dt, 

I=—F{X) j nd{r,a)d(l>. 

Mathis, Rumpl, & Nordsieck (1977, hereafter MRN) found exctinction of starlight through many lines 
of sight to be well reproduced if the distribution of dust grains in the ISM (for a given range of a) is given 
by 

dugr = Cniia~^'^da (6) 

where ngr{a) is the number density of grains with radius < a, C is a normalization constant, and na is 
the number density of H nuclei. Since dugr/da and nd{r,a) are both the number density of grains between 
(a, a+da), equation (6) can be rewriten nd{r, a) = nH{r)f{a) where the exact form of /(a) will be determined 
in §2.4. 

Light echoes arc observed as extended flux distributions, therefore we measure the surface brightness 
Bsc of scattered light. Noting that Fsc = / BscpApdcj) and Fsc{X,t) = / Fsc{X,t,(j))d(j), equation (4) now 
becomes 

Bsc{X,t,cl>)=F{X)7i^{r)^—^^S{X,p) (7) 
where the integrated scattering function 

S{X,n)= j Qsc{\a)ag^{n,X,a)f{a)da. (8) 

Since the observed flux F{X) and area element pAp diminish as D~'^, equation (7) for surface brightness is 
distance independent, as expected. If the light echo is close to the line of sight to O, the echo should suffer 
roughly the same extinction as the observed flux F{X), mitigated of course by the potentially higher dust 
density directly surrounding the echo location. 
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2.2. Light Echo Width 

Since a light pulse has a finite duration r, the observed light echo will have a width dpg resulting from 
the convolution of the dust thickness A2: and the pulse width 

cz + c?t 

dpt = T. 9 

P 

This is shown schcimatically in Figure 2. Xu at al. (1995) approximate the echo paraboloids in the neighbor- 
hood of i^z by lines of slope k = dp/dz, and find the convolution of two gaussians of FWHM oc A^r and drt 
yields 

{kl2dp}f 

2dpl 2dp1 (fcV2dp2) + (i/2A^2) 

which simplifies to 



,2 



A.= vMEM. (10) 



Using the light echo equation [eq. (1)] one finds 



- = -,anddp, = ^- + -)r 



with which equation (10), solved for Ap = dpe, becomes 



Since z and r = \J p^ ^ z^ are both functions of (p, f), the geometric factors in equation (7) 

and p, can be written in terms of observables (p, r) with A2; the only unknown. This allows an echo to be 
specified in terms of the dust structure (Az) rather than its observed width (Ap). The geometric factors G 
and Q arc shown in figure 3. In discussing the geometry it is often easier to use units of years and light-years, 
since c — I. 



2.3. Time-Integrated Flux 

To compute the surface brightess of an echo, one must know the time-integrated monochromatic flux 
F{X). If a homogenous time-series of spectra exist for a given object, one simply integrates them over time, 
provided they have sufficient spectral coverage for the wavelength range in question. For this paper, I wish 
to explore the observability of echoes around a variety of objects, and over a large wavelength range (1000- 
11000 A). Time -series of spectra are available for many (but not all) objects of interest, and the wavelength 
coverage is highly inhomogenous, thus such an approach for determining F(X) is not ideal. Rather, I propose 
a method by which a single spectrum can be used to approximate the time-integrated flux. This has the 
added bonus that it defines the value of r [eq. (12)] in a consistent manner. 

Although a light pulse has finite duration, its flux varies substantially during r. Furthermore, deter- 
mination of the exact value of r is ambiguous - some commonly-used values are the time for flux to vary 
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by two or three magnitudes, or the "effective width" J F{t)dt/ F{tmax)- Since fainter fluxes contribute less 
to light echo signal, a compromise must be reached between inclusion of a reasonable portion of the light 
curve (i.e. the wings), and exclusion of faint fltixes whose contributions to an echo are negligeable compared 
to the peak flux. Sample light curves of composite Type II-P SNe (Doggett & Branch 1985) and the Mira 
star X Cyg arc shown in Figure 4. These curves can be approximated in the neighborhood of their maxima 
by parabolae in magnitude space. A parabola with vertex (0, log/o) and passing through (i',log/i) can be 
written 

log/(t)= log/o -^L (13) 
where L = (log/o — log/i). The area A under the curve between ±t is found by integrating /(<): 

A{t) = 4^erf(tx/LhaO/i'). (14) 
VLm 10 



The error fimction erf approaches unity if the operand is > 1.4, that is if t > tc = \At' j^jL InlO. Let 
tn be the time required for the flux to drop by n magnitudes, or equivalently, for log[/(t)] to drop by 0.4n. 
Applying this constraint to equation (13), one finds i„ = f'-\/0.4n/L. Equation (14) is roughly constant 
for t„ > or n > 2.1. Thus, a reasonable value of r is 2*2, i.e. the time required for the flux to rise then 
diminish by 2 magnitudes. At r = 2*2, equation (14) reduces to 

By choosing to construct the parabola with t' = t2, VL = \/0^ and A(*2) = 1.25/o*2. With this, an 
approximation for the time-integrated monochromatic flux F{X) in equation (7) can be expressed by F{X) = 
1.25F(A, t„i,ax)t2-, where tmax is the time of maximum light. 

How does this prescription for measuring r compare with other methods? To address this, I have 
calculated a number of diagnostic widths of the mean B light curves for type la, II-P and II-L supernovae 
from Doggett & Branch (1985), which the authors provide in an interpolated tabular form. The results are 

given in Table 1. The left columns give the values of r for the five methods noted below. The right columns 
show how much of the total flux of the light curve is included in this measurement of t. The methods are: (1) 
the effective width reported by Sparks (1994); (2) the effective width I calculate integrated from maximum 
light; (3) the effective width I calculate by integrating the whole light curve; (4) the time for the light curve 
to actually rise and fall by 2 mags from maximum light; (5) twice the time for the light curve to fall 2 mags 
from maximum light. There is a large dispersion among the different methods, demonstrating the ambiguity 
of this measurement. The effective width measured from maximum light is roughly half the value to rise/fall 
by 2 mags, showing that despite the dispersion among methods, they are qualitatively consistent. Aside 
from the type II-L curve, whose integrated flux saturates very quickly, the effective-width measurements 
(methods 1 3) exclude roughly one-quarter of the energy released in the outburst. 

To some extent, it is up to the individual to select which prescription for r is preferred. In the context 
of this paper, the *2 formulation presented above gives an estimate not only of the width of the light pulse, 
which is needed in G [eq. (12)], but also of the time-integrated flux F{\) that enters equation (7). For 
the purposes of exploring a general parameter space, I will use equation (15), and caution the reader that 
integrated fluxes may be in error by 10-25%. In practice, the spectral energy distribution of a variable object 
changes over r, which this prescription for F{\) does not take into account. Generating an empirical time- 
integrated spectrum is the obvious method of choice, provided (as noted above) one has sufBcient spectral 
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and temporal coverage. In §4, I discuss the error associated with this approximation apphed to SN 1993J, 
and present an easily-implemented method of correction. 

2.4. Dust Properties 

In their classic analysis, MRN found that a mixture of carbonaceous and silicate grains following a 
density distribution given by equation (6) adequately modeled observed starlight extinction. Draine & Lee 
(1984) and Laor & Draine (1993) calculated the dielectric functions for graphite and "astronomical silicate," 
from which they calculated the dust properties (QsC) QabSj and g) using the dipole approximation (valid 
for a/A < 0.1) and Mie theory (van de Hulst 1957). Qabs is the absorption efficiency, relevant when 
considering IR emission or optical extinction. Weingartner & Draine (2001, hereafter WDOl) have corrected 
an anomolous feature in the silicate dielectric function, while Li & Draine (2001) have included the effects 
of small PAH molecules in carbonaceous grains. Values of Qsc(A, a) and g{X, a) are shown^ in Figure 5, for 
both carbonaceous and silicate dust over a variety of grain radii and wavelengths. 

WDOl have constructed a new dust-grain density distribution /(a) fit to observed interstellar extinction. 
I refer the reader to their paper for the actual functional form, but note that it contains a smooth cutoff for 
large and small grain sizes and allows the slope d\nngr/d\na to vary, unlike the fixed slope of equation (6). 
I have adopted their preferred Milky Way dust model for {case A, Ry = 3.1, and 6c = 6 x 10~^}, shown in 
Figure 6, which works well at fitting both optical extinction and infrared dust emission (Li & Draine 2001) in 
diffuse clouds. For comparison, more reddened sightlines are well fit by the {case A, Ry = 5.5, &c = 3x 10~^} 
density model, also shown in Figure 6. The reddening ratio Ry = Ay / E[B — V) (Cardelli, Clayton & Mathis 
1989), and 6c is the total C abundance per H nucleus. I caution that these particular density distributions are 
not unique, and WDOl discuss a range of parameter space that fit the observed galactic optical extinction 
(roughly) equally well. Furthermore, extragalactic dust densities need not follow that of our galaxy; the 
authors find a different dust distribution for the LMC extinction, for example. Uncertainty in the density 
distribution is to be expected, since e.g. (as noted by WDOl) the dielectric functions are not precise, dust 
grains have poorly-understood surface monolayers, and grains are nonspheroidal. As such, predictions based 
on the adopted interstellar dust grain distribution should be taken with a proverbial salt grain. 

2.5. Exposure Time Calculation 

Since the aim of this paper is to discuss light echo observability, I briefly review the calculation of 
expsoure times for a generic telescope, instrument, observing conditions and input spectrum. The signal-to- 
noise {S/N) ratio for an exposure time t is given by the CCD equation (Howell 2000) 

S/N = Coi,jt/[{Cobj + C,ky + Cbg)t + nocit + tno) + nlj, + G^a}]^^ (16) 

where tno is the read-out time of the detector; G is the CCD gain [e~ DN~^]; 1 DN is one "data number," 
or count (also known as ADU); noc is the dark current [e~ s~^ pix~^]; nun is the read noise [e~]; cr/ is 
the digitization noise of the digital-to-analog converter [DN] and is ^ 0.3 (Howell 2000); and C indicates 
the count rate [e~ pix~^], with Cobj for the object observed, Csky is that of the sky, and Cbg is the 



Tabulated values are available at http://www.astro.princeton.edu/~draine/dust/dust.diel.html 
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background, e.g. nebulosity or galaxy surface brightness. Dark current, read time, gain and read noise are 
properties of the detector. For a desired S/N it suffices to solve the quadratic in equation (16) for t. 

The count rates are computed from the sensitivity of the system (as clearly explained in Leitherer et 
al. 2001), which gives the counts per unit wavelength resulting from unit incident flux. The sensitivity is 
calculated as (the effective area of the telescope at A) divided by (the energy per photon at A). Let us define 
Tx(A) to be the total throughput of the observational configuration X (e.g. telescope, filter, detector) at 
wavelength A, then sensitivity is Tx{^)AX/hc where A is the area of the primary mirror, and he is Planck's 
constant times the speed of light. It follows that for an object with surface brightness B{X), the count rate 



where mx is the platescale [arcsec pix"-'^] and Ax is the extinction in filter X. For a 2.4-m telescope (typical 
of a moderate-sized ground-based observatory or HST) A/hc ~ 2.3 x 10^^ cm^ erg^^ A^^. Throughputs for 
ACS, STIS and Johnson filters, including HST optics, are taken from synphot. Although a ground-based 
telescope will not have the same throughput as HST, this approximation should be sufficient. For HST, 
detector characteristics and "typical" sky brightnesses (my = 22.7 mag arcsec""'^, including zodiacal and 
geocoronal light) are taken from the Instrument Handbooks (Pavlovsky et al 2002; Leitherer et al. 2001). 
For a generic 2.4-m ground-based telescope, I adopt sky brightness estimates for 3 days after new moon from 
Walker (1987), and use characteristics typical of a SITE 2048 x 2048-pixel thinned and backside-illuminated 
CCD with 0.24/im pixels, with platescale computed at f/15. 

For ground-based observations, I consider Johnson and Cousins filters between U and /. STIS has 
three cameras of interest: the optical long-pass F28x50LP filter and the two MAMA detectors (near and far 
UV). For ACS, I consider the optical Wide Field Camera (WFC), the optical and near-UV High- Resolution 
Camera (HRC), and the UV Solar-Blind Camera (SBC). Characteristics for all filters and detectors used in 
this work are listed in Tables 2 and 3, respectively. 



As noted in §2.1, an echo and source should suffer the same extinction, provided their angular separation 
is small. Using the observed (i.e. not extinction corrected) source spectrum for -F(A), extinction need not 
enter equations (7) or (17). Since in this work I will often use template spectra (§3.2), extinction must be 
included in surface-brightness calculations. Furthermore, it is important to know how much extinction is 
caused by the echoing dust, in order to rule out unreasonable regions of the parameter space (i.e. density or 
composition). 

The power per unit area of an electromagnetic wave traveling a distance L through dust of density na 
will be attenuated by a factor cxp (— n^iCEXT) (Draine & Lee 1984). The factor Cext is the extinction 
scattering cross section, defined as the sum of the scattering and absorption cross sections Csc and Cabs- 
Since the cross section is related to efficiency hy Q = C/ag, the extinction (in magnitudes) can by calculated 




(17) 



2.6. Extinction 



by 




(18) 



where A^h = tihL is the column density of hydrogen nuclei along the optical path. 



I have integrated equation (18) using the dust density, scattering and absorption properties discussed in 
§2.4. Integration was performed using the extended Simpson's method (Press et al. 1992), with values for Q 
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and g interpolated logarithmically in A and a. Figure 7 shows the resulting extinction for Ry = 3.1 and 5.5 for 
Galactic dust, as well as the individual contributions from silicate and carbonaceous dust. For comparison, 
I also show the extinction curves from Cardelli et al. (1989) in each panel. As these are reported as Ax/ Ay, 
I have normalized them to the computed total extinction at 5495 A (the value at which Ax /Ay = 1). The 
fit is very reasonable for both redcnning parameters. From these curves, the value of Ax in equation (17) 
can be calculated, given either an estimate for Ay or A^h- 

2.7. Non-Planar Dust Distributions 

Figure 8 shows schematic situations typical of a CSE around a mass-losing star. Mass-loss produces 

an (gas) density profile, shown here as a grcyscale gradient. The inner and outer radii are chosen such 
that Rin designates the radius at which dust condenses in the wind, and Rout is a matter boundary, such as 
when the wind density reaches that of the ISM. As discussed in §2.2, one will observe an echo from all the 
material located between the echo parabolae marking the beginning {t{) and end (t2) of the light pulse or 
the star's high-flux state. An echo in a CSE differs from the simple case explored in the previous subsections 
since density now varies as a function of scattering position. Some variable stars have periods larger than 
the light-crossing time of the envelope, thus the envelope material swept up between the echoes at ti and t2 
will be substantial, as shown in panel (a). In this "thick echo" regime, one may consider the light echoing 
from a single column of width Ap of this envelope. Although the overall geometry is not planar, one can 
apply the planar approximation by integrating thin planar contributions along this column. If the period 
of variation is smaller than the envelope size, one may consider the "thin echo" situation depicted in panel 
(b), where a single sight line looks through many nested parabolae (for clarity I have only drawn 12 echoes, 
although 39 fit into this geometry). The contributions from all such echoes must be included. 

If the envelope has high optical depth, the single scattering approximation is no longer valid. For the 
purposes of this work, single scattering will be sufficient, since Ay < 1 in most dust geometries (§3.3.4). A 
more thorough analysis of multiple scatterings should be investigated for cases when Ay > 0.3, however this 
first approximation will serve as a guide to whether such echoes will be observable at all. 

3. THE LIGHT ECHO MODEL 
3.1. Dust Scattering 

Computation of the integrated scattering function ^(A, ^) is a straightforward integral over grain sizes, 
using equations (3), (8) and the dust-density models from WDOl. Selected results are plotted in Figure 
9, showing the dependency of 5(A,/x) on density distribution model, grain size range, composition, and 
scattering angle. Columns show a given scattering angle 6 with varying grain-size range, and rows show a 
fixed size range for varying 6. Changes along columns are due strictly to the limits of integration for S*, 
while changes along rows are the result of variations in the scattering phase function The three size-range 
distributions are denoted "S" (small, a = 5 x 10~^ - 0.01/xm), "M" (mid-sized, a = 0.01 - O.l/^m) and 
"G" (Galactic range, a = 5 x 10~^ — l.O/im). "Galactic" dust is the full density distribution adopted by 
Weingartner & Draine (2001). 

Scattering for grains with a < 0.1A/27r occurs within the Rayleigh limit, for which S oc A"**. This 
behavior is clearly seen in the top row in Figure 9 and in the second row, for A > 5000A. For the smallest 
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silicate grains, S oc A~^'^, while 5* cx A~^'^ for carbonaceous dust with A > 3000A. This "super- Rayleighan" 
behavior reflects the dielectric function and resonance effects (B. Draine 2002, private communication). As 
the maximum grain-size increases, scattering efficiency increasingly deviates from A~^ behavior, but UV 
scattering is generally more efficient than optical. 

Of particular interest are the scattering peaks seen in the UV (1000-2500A). The presence of the strong 
1500A carbonaceous peak/dip in the bottom row (type G dust) demonstrates the tendency of large dust 
to predominantly forward scatter near this wavelength. This conclusion is borne out by observing how the 
scattering function changes with n — cos 6. Regardless of the adopted form for the scattering phase function, 
$ should be larger (smaller) for scattering toward (away from) the observer. For equation (3), as /j. — > ±1, 
$ ^ (1 ± .g)/(l T .9)^7 thus as the forward scattering efficiency g approaches 100%, 3> grows large for 6 — 
and small for = 180°. The degree of forward scatter g > 0.8 for larger grains (a > O.lAtm) at 1500A, but 
g < 0.04 for grains smaller than 0.01/um. Large grains tend to forward-scatter UV light, while small grains 
scatter more isotropically. Thus we expect a UV bump near 1500A for ^ ~ large-grain scattering, and a 
dip for 6 — > 180°, as observed in the bottom row of Figure 9. 

Maxima in S near lOOOA and 2175A in carbonaceous grains arc clear when the large-grain A = I5OOA 
resonance is not present, i.e. for the smallest grain sizes and for — * 180°. These two features demonstrate 
the high scattering efficiency Qsc of dust at short wavelengths. The 2175A bump is a known interstellar 
scattering and extinction band (Stecher 1965; Stecher & Donn 1965), which WDOl attribute to PAHs. The 
loooA maximum is an artifact of the wavelength range plotted, and S continues to increase for shorter 
wavelengths (there is another bump around 800 A). Large silicate grains also show a bump around 2000 A, 
the origin of which is unclear. 

Figure 9 suggests that a great deal about dust can be learned from carefully-planned observations of 
light echoes. The actual light echo flux will be a function of the total dust density, grain size distribu- 
tion, and the relative population of carbonaceous-to-silicate grains, geometry and light-pulse characteristics 
notwithstanding. Since we generally know the latter two, echo colors and magnitudes should constrain the 
dust properties. Assume for now that the density distribution function is known, and only the grain-size 
limits, density and composition are not. The observed echo spectrum will be a product of the dust scattering 
S{X,fj.) and the pulse spectrum F{X) (§3.2), from which object-specific echo models can be created. 

Echo geometry fixes indicating on which column in Figure 9 the total dust scattering efficiency must 
lie, while the dust model can be chosen from Ry, as determined from stellar spectrophotometry. Optical 
colors alone may distinguish the largest from the smallest dust grains (since smaller grains at longer wave- 
lengths follow the Rayleigh scattering relationship) , however since the dust number density (expressed here in 
terms of hydrogen density) will likely be unknown, brightnesses themselves will not conclusively distinguish 
between various particle-size scenarios. Fortunately, blue to UV fiuxes are a strong measure of dust size, 
and UV-to-optical colors should lift the small-particle degeneracy. Furthermore, aptly chosen wavebands 
can distinguish emission from the various scattering maxima, which are sensitive to dust composition. For 
example, the STIS MAMA far and near UV filters have average wavelengths and FWHM of (138lA,324A) 
and (2310A,1237A), respectively, with only mild overlap at 1500A. The question remains whether echoes 
at these wavelengths are easily observable. Finally, it is important to note that spectroscopy, even in low 
resolution but covering a wide wavelength range, is likely the best observational strategy for studying dust 
properties through light echoes. 
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3.2. Variable Objects: Selection and Spectra 

To predict the brightness of echoes from different sources, it is first necessary to select candidate objects 
to study. The primary selection criterion is that the object vary in brightness by a significant amount. 
Without variability, any scattering material will appear as a reflection nebula rather than as echoes, and no 
3-D information can be gleaned. A second criterion is that the duration of variability is short compared to 
the light-travel time across at least one dimension of the scattering medium. A burst of duration 100 days 
can not give any spatial information about a dust shell of radius 1 light day, but can still map the dust 
in a thin but extended dust sheet since the sheet's dimensions are wider on the plane of the sky than the 
pulse width. I list in Table 4 classes of variable objects that change in brightness by AV > 2 mag. For 
easy comparison, the data are also shown in graphical form in Figure 10. Four broad categories are present: 
supernovae, cataclysmic variables (CVs), pulsating giants and eruptive stars. 

To model echo brightness, and to discuss the relative brightness between different wavebands, it is 
necessary to provide a realistic representation of the object's emission spectrum in the factor F{X) in equation 
(7). Selected representational spectra for the first three categories are plotted in Figure 11, and discussed 
below. Each spectrum has been normalized to a distance of 1 pc. This slightly-unconventional normalization 

has been chosen to balance the conversion from physical area {pAp)~^ in G (§2.2) to angular area. Since 
these will represent the flux at maximum brightness, I will refer to these as F{X.trnax)- 

Supernovae: Immediately following core collapse, type II supernova shock breakout is brightest at 
shortest wavelengths. The short-lived but intense pulse of UV emission will photoionize any surrounding 
circumstellar material, and likely evaporate the smaller, nearest (d <C 1 pc) dust grains (see Emmering & 
Chevalier 1989), however larger and more distant grains should survive. As shown in §3.1, dust scattering 
is most efficient in the UV for most echo geometries, making dust echoes of the UV pulse an attractive 
observational target. The later- time spectrum is dominated by optical emission from cooling lines and 
radioactive decay, and provides the longer-duration optical pulse that is reflected in optical light echoes. To 
model echoes from SNe, I use synthetic spectra of SN 1987A from day 1 after core collapse (when the SN 
was UV bright; Figure 11a), and day 58 after core collapse (nearing the maximum optical light; Figure lib), 
generated with the Phoenix stellar and planetary atmosphere code (E. Baron 2002, private communication). 
Both spectra were normalized using a LMC distance of 50 kpc. A Phoenix synthetic spectrum of a Type la 
event 25 days after explosion (not shown) has stronger flux than SN 1987A (a Type II SN) at day 58, and 
while the microscopic features are different, the overall spectrum is similar. I will therefore use the Type II 
spectrum for echo calculations, with a caution to the reader that type la events are 1-2 mag brighter. 

CVs: Within this category are classical, recurrent and dwarf novae, and symbiotic stars, together 
spanning up to 10 mag between the brightest and faintest outbursts, with such mechanisms as thermonuclear 
runaway on an accreting white dwarf to accretion disk instabilities. To model echoes from classical novae, 
I use synthetic spectra of a nova outburst during the fireball (Figure lie) and constant luminosity (Figure 
lid) phases, again generated with the Phoenix code (P. Hauschildt 2002, private communication). Both have 
been normalized to Nova LMC 1991 (Schwarz et al. 2001) assuming an LMC distance of 50 kpc. Recurrent 
novae (including U Sco and RS Oph) should be adequately modeled as classical novae with slightly less flux. 
Dwarf novae (DNe) come in a variety of flavors: U Gem, Z Cam, SS Cyg and SU Uma DNe have amplitudes 
of 2-6 magnitudes during normal outbursts, while TOADs - SU Uma and WZ Sge categories - can have 
outburst amplitudes of 6-10 mag (Howell, Szkody, & Cannizzo 1995). As DNe are abundant, quasi-periodic, 
and can be bright, they may provide a generous target list of echo candidates. 

DNe occur from instabilities in binary-star accretion disks. The theoretical spectrum of an infinitely- 
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large, steady-state accretion disk is given by F\ oc X~'^^^ (Lynden-Bell 1969), however since the disk has a 
discrete inner edge (the accreting white dwarf) this law is more reasonably written Fx oc A~^-^ (J. Patterson 
2002, private communication). A normalized (to Ipc) spectrum is easily generated from this law using SS 
Cyg: F(A) = lAX''^ '^ where A is in A (Polidan & Holberg 1984; Harrison et al. 1999). 

Pulsating Giants and Supergiants: Pulsating giant stars include (but are not limited to) long-period 
variables (LPVs) such as post-ABG thermal-pulsating carbon stars and Miras, G-K type supergiant RV 
Tauri stars, and F-K type Ccphcids. Figure lie shows the spectrum of the M7III star SW Vir, available in 
the Bruzual-Persson-Gunn-Stryker (BPGS) spectrophotometry atlas within synphot. Itself a semi-irregular 
variable, SW Vir has a spectral type similar to the prototypical Mira o Get. This spectrum has been 
normalized using the star's Hipparcos parallax of 7 mas (Ferryman et al. 1997). As expected of a cool giant, 
the energy distribution is predominantly red. 

A Ccphcid has a type Fib spectrum during maximum brightness. Many such stars exist in the Jacoby- 
Hunter-Christian (JHG) spectrophotometric atlas within synphot, but only between 3500-7400A. To generate 
a Cepheid-like spectrum with the full wavelength coverage needed, I scaled the FOIV spectrum of ^ Ser (from 
the BPGS catalog) to the F6Ib spectrum of HD 8992 (from the JHG catalog), having first normalized HD 
8992 to a Hipparcos distance of 0.82 mas, and corrected ^ Ser for an extinction of Ay = 0.015 (listed in the 
file header). Figure llf shows the resulting spectrum. 

Eruptive Stars: Many stars in this class do not satisfy the required criterea, such as R Goronae Borealis, 
which fades rather than brightens; UV Ceti stars, which have pulsation durations of seconds to minutes; or 
Wolf Rayct stars, which only vary by AV -^0.1 mag. However, the S Doradus class is a high-luminosity 
blue star, typically surrounded by diffuse nebulosity or expanding envelopes, rj Gar is a famous member of 
this class, but with a period of over 5 years. I note this class for completeness, and although echoes might 
be observable in the extended envelopes, I will not make predictions for these stars. 

3.3. Observability of Light Echoes — Constraints 

Prior to integrating equation (7), we must specify the echo geometries and dust densities for each 
candidate object. Once specified, this integral yields the source count rate to be inserted into the GGD 
equation [eq. (16)], along with the appropriate values from Tables 2 and 3 and an estimate of the background 
Cbg ■ An echo is deemed "observable" if it can be detected at a given signal-to-noise level in a "reasonable" 
exposure time. In the following sections, I discuss estimates for exposure time, dust density and geometry, 
extinction, and sources of background flux. 

3.3.1. Exposure Time 

Since light echoes are transient objects, they are effectively detected by PSF-matched difference imaging, 
in which one matches the PSFs between two images using a variety of techniques prior to subtracting them. 
This offers a minimally-invasive method for removing all sources of constant flux, with only sources of variable 
brightness remaining. The differencing technique of Tomaney & Crotts (1996) has been shown by Sugerman 
& Grotts (2002, hereafter SG02) to be quite effective at detecting light echoes in HST imaging, as well as 
in e.g. the monitoring of hot spots in SNR 1987A (Lawrence et al. 2000; Sugerman et al. 2002), detection 
of variable stars in globular clusters (Uglesich et al. 1999), and the M31 microlensing project (Grotts 1992; 
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Uglesich 2002). 

Direct photometric techniques (such as aperture photometry or even PSF-fitting and removal) can only 
remove resolved point sources, but have difficulty identifying faint surface-brightness features within the 
fluctuations of an unresolved background component. In differencing two images, these unresolved sources 
are removed, leaving (in principle) only statistical noise. Even though this noise increases by \/2, the 
detection tlircshhold (i.e. required S/N) for an echo decreases notably in a difference image, since we are 
looking for coherent signal over a large, and now statistical-noise dominated, area. A well-resolved echo need 
only be detected at S/N = 2 — 3 (per pixel) to appear unambiguously above the difference-image background. 
I will adopt a per-pixel S/N = 3 as the required detection threshold. 

The assessment of a "reasonable" exposure time is subjective, but not without practical constraints. 

A reasonable series of HST observations would image in 3-5 filters within perhaps double the number of 
orbits. With 50-70 minutes per orbit, a 30 90 minute exposure seems a reasonable maximal exposure time. 
Ground-based observations are much more flexible, and I will consider a six-hour (total) exposure time in 
any given filter as a representational (but by no means fixed) upper limit. 



3.3.2. Gas and Dust Density 

For all model predictions, I will assume dust follows the type G (or Galactic) dust distribution (§3.1). To 
calculate the dust density, it is first necessary to specify nn, the density of hydrogen nuclei, which normalizes 
dugr/da to give the dust density within interstellar gas. The use of this normalization factor depends on 
the assumption of a constant gas-to-dust ratio in the ISM. however this same ratio may be unreasonable 
for circumstellar environments. The dust distribution integrates to yield a gas-to-dust ratio of 124, and 
total number densities (per material) of nc = 6.0 x 10~''nH and risi = 3.5 x 10~^nH for carbonaceous and 
silicate dust (respectively) with Ry = 3.1, while the Ry = 5.5 distribution gives nc = 3.0 x 10~^nH and 
nsi = 4.5 X 10~^°nH. For a given dust composition {nd = one + ^'T'Si), the above densities can be used to 
rescale nn to the appropriate dust content within any environment. 

Comparison of these dust densities with previously-reported values for the ISM will be problematic, 
since the WDOl density profile favors a large number of extremely small grains. Small particles are inefficient 
scatterers at optical wavelengths: Qsc = 10"^*^ — 10~^ for particles with a = 3 x 10~^ — W~^^m, while 
Qsc > 10~^ for a > 5 X 10~^/im (Fig. 5). These grains contribute little to optical scattering, and have 
been ignored in previous density distributions. Their effect on the total density is easily understood by 
considering the MRN distribution in equation (6), which integrates to yield 

This behavior strongly biases the density by their presence, yielding values much larger than generally quoted. 

For a meaningful comparison with previous work, I have integrated the WDOl distribiition over the size limits 
given by MRN: 0.005 < ac < 0.25 and 0.01 < asi < 0.25. With these limits, the total number densities for 
Rv = 3.1 are nc = 7.6 x 10~^^nH and nsi = 2.3 x 10~^^nH, while for Ry = 5.5, nc = 4.1 x 10~^^nH and 
nsi = 2.8 X lO-^^nn. 

Interstellar Dust: Following the multi-stage ISM model of McKee & Ostriker (1977), the majority of 

interstellar gas will be in either the hot (ud > 10~^) or warm (n^ ^ 0.1 — 1) states. Were we to regularly 
observe echoes from this ambient-density medium, one would expect to see a diffuse halo around echoing 
objects, rather than discrete arcs and rings. This suggests that we detect most light echoes from only the 
denser regions of the ISM, however molecular clouds (nn ^ 1000, r ^ 50pc) have sufficiently high extinctions 
{Ay > 5) that echoes should not be observed. Rather, echoes likely originate from HI clouds with nu ~ 10, 
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which have sheet-like features and extinctions Ay < 1 (R. Chevaher 2003, Private Communication). 

Mass-losing giants: Post-AGB stars have mass-loss rates of 10~* to 10~^ Mq yr^^, yielding circumstcUar 
density profiles = Kr~^ (where if is a constant of proportionality) , and with dust condensing within the 
outflow as near as 8i?*, where i?« is the stellar radius ~ 3 x 10^"^ cm (Habing 1996). Justtanont, Skinner, & 
Tielens (1994) find gas-to-dust ratios of ~ 200 for oxygen-rich red- giant winds, which implies that a factor of 
0.6nH can be used to predict the dust density. Vassiliadis & Wood (1993) find for Miras a mass-loss scaling 
relation logM = -11.4 + 0.0123P(days), yielding M ^ 10"^ Mq yr^^ for a 200 day Mira. The equation of 
continuity predicts that the constant if = 3 x IQ'^^M /voo where M is in Mq yr~^ and v^o is the gas terminal 
velocity in km s~^, typically 5-25 km s~^ for AGB stars, and 30 km s~^ for cool giants (Lamers & Cassinelli 
1999). As such, ff ^ 3 x 10^"^ cm"-'^ for 200-day Miras. Ccpheids have mass-loss rates between 10~^° to 
10~^ M© yr~^ (Deasy 1988). While most cepheids shed mass at the slower rate, high mass-loss examples 
are known, such as RS Pup, whose mass-loss is estimated by Deasy (1988) to be ^ 10~^. At this rate, the 
density law will have K = IQ^^ cm-\ while at M = 10"^, K = 10^^ cm'^ 

CVs: Dust can condense in a stellar wind once the equilibrium temperature of the gas is less than 
^ 1500K, which occurs (Lamers & Cassinelli 1999) at a radius 



where p ~ 1 for silicate and p 2 for carbonaceous dust. For a dust condensation temperature ~ lOOOK, 
p — 1.5, a white dwarf/accretion disk effective temperature of (4— 10) x 10'' K. and an inner radius 7?* — 10^ 
cm, this yields Tj, = (1 — 16) x 10^^ cm. UV resonance-line profiles from DNe imply mass loss rates of 
M = 10-9 Mq yr-i (Vitello & Shlosman 1993) with terminal velocities in the range of 3000-5000 km s"^ 
(Warner 1995). For the smaller- velocity limit, the continuity equation gives K = 10^^ cm"^. 



The geometric parameters entering G [eq. (12)] are listed in Table 5 and discussed below. 

Mass- Losing Giant Stars: Post-AGB and Mira stars (and realistically, most stars late in their lifetimes, 
such as rj Car, SN 1987A) have complex stellar outfiows, often bipolar or asymmetric, and are of critical 
interest as progenitors of an increasingly complicated taxonomy of planetary nebulae (as demonstrated in 
Kastner, Soker & Rappaport 2000). Although the sc;attc;ring geometry is not planar in dusty envelopes (as in 
Figure 1), as discussed in §2.7, we can still apply the planar approximation. As shown in Figure 8, the echo 
will occur between two parabolae at ti and t2, separated in time by r. The choice of geometric parameters 
representing typical echoes is somewhat arbitrary. Below, I discuss the geometries I will test. 

The prototypical Mira x Cyg has an inferred CO shell extending to 3 x 10^^ cm (Jorissen & Knapp 
1998), while the post- AGB carbon star V Hya shows evidence of bipolar outflows at radii of ~ 1.5 x 10^^ 
cm (Sugerman, Sahai & Hinkle 1996). As the pulsation periods are large compared to this size, I use the 
thick-echo approximation (§2.7). As suggested by the aforementioned stars, one radius to test is 3 x 10^^ 
cm, as well as 10^^ cm. For both cases, I choose ti so an echo passes through r with 6 = 30°. I take the 
period to be 200 days, r=50 days, and M = 10~^ Mq yr~^. This is shown schematically in Figure 8a. 

For Ccpheids, I use an envelope geometry suggested by RS Pup. Variations in light intensity in the 
surrouding nebula were interpreted by Havlen (1972) as arising from shells of reflecting material spaced 




(19) 



3.3.3. Dust geometry 
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between inner and outer radii of roue hly 5 X 10" and 2.3 x 10^^ cm. This inner radius is roughly ^ly, much 
longer than any Cepheid pulsation rate, thus the thin-echo approximation applies. I consider the innermost 
echo to be located within a shell with radii given above, with ti such that the inner echo passes through 
Rout a.t p = Rin- This is shown to scale in Figure 8b. I test two cases: that for the high mass- loss rate of 
RS Pup, and a more typical Cepheid mass-loss rate of 10~^ Mq yr~^. 

Classical Novae: The first nova outburst will eject 10~^ — 10~^ Mq of material at hundreds to thousands 
of km s~^ (Bode & Evans 1989), expanding into the ISM in a manner analogous to that of a SN, the physics of 
which is well studied (e.g. Chevalier 1974). Indeed, this is believed to be the case for Nova Per 1901 (Seaquist 
et al. 1989). Light echoes from this event may highlight the CSE into which the ejecta will travel, thereby 
allowing us to better understand (in the future) how the ensuing nebula formed, as well as the mass outflow 
characteristics of the star pre-nova (i.e. in its quiescent binary period). Recurrent novae may echo off their 
swept-up ejecta shells, providing direct testing of predictions for these lower-energy ejecta. Furthermore, 
such echoes will constrain the recurrence time of novae, a quantity which remains only speculative. Both 
types of bright novae may illuminate undisturbed ISM material as well. 

Nova Per 1901 was a bright {AV ^ 13) classical nova, and notably it is the flrst object aroimd which 
light echoes were detected. Seaquist et al. (1989) have analyzed the archival imaging of GK Per following 
its 1901 outburst, concluding that the echoes analyzed on day 210 by Couderc (1939) were from a plane of 
dust with r = 1.4 pc at nearest approach, and with thickness Az ~ 0.15 pc. It is unclear if this was true 

interstellar material, or a previously swcpt-up shell of circumstellar material. I will consider a sheet of dust 
with these characteristics, as well as dust with (somewhat arbitrary) thicknesses of Az = 5 ly at z =50, and 
Az = 25 ly at ^ = 250 ly. 

Dwarf novae: Being much less energetic in their outbursts than classical novae, DNe may not evacuate 
their circumstellar material in the same way, and hence may illuminate their quiescent stellar mass loss. GK 
Per currently has DN outbursts, the brighter of which might echo off the swept up nova shell. Of particular 
interest is the geometry of outflows near the white dwarf. While most accretion systems drive narrowly- 
collimated jets (e.g. T Tauri stars. X-ray binaries, quasars), non-magnetic CVs show evidence for mass loss 
with large inferred opening angles (Knigge & Drew 1997), suggesting these are accretion systems in which 
jets are absent (Knigge & Livio 1998). I consider a TOAD nova echoing from circumstellar (r = 0.025 ly) 
and interstellar (r = 25 ly) dust, as well as a dust sheet similar to that just discussed from N Per 1901. As 
the thickness of dust within the DN wind is uncertain (due in part to the unknown orientation of the outflow 
axis) I will arbitrarily take Az to be a light day. 

Supernovae: Supcrnovae offer the means to ilhmiinatc the most distant material, thereby probing the 
largest portions of the ISM. Unfortunately, the observability of resolved echoes around SNc is limited by 
the low frequency of nearby events. As in the previously-discussed objects, echoes from circumstellar and 
interstellar dust offer invaluable information on the progenitor's evolution and the surrounding ISM. In 
particular, echoes can probe multiple stages of progenitor mass loss as imprinted in "snow-plow" contact- 
discontinuities with the ISM (e.g., Crotts & Kunkel 1991, for SN 1987A). I study six possible dust geometries, 
two each from circumstellar, contact discontinuity, and interstellar regions, using geometric parameters 
similar to echoes observed around SN 1987A (Crotts et al. 1995; Xu et al. 1995). 
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3.3.4- Extinction 

To account for interstellar extinction from the light echo to Earth, I use the extinction curve in Figure 

7, which relates A\ in terms of Ay- However, the choice of this parameter depends on the particular line of 
sight and the distance of the echo from Earth. Rather than consider the myriad sightlines of objects, I adopt 
Ay = 0.25, with the caveat that predictions for specific objects may require different extinction values. 

Mass-Losing Giants: Consider the trajectory of photons from the central star to the point of scattering 
and subsequently out of the CSE. In the single-scattering limit, light is extincted along both paths. An 
estimate of the extinction can calculated using the extinction curve in Figure 7, from which Ay = 3 x 10~^^iVH 
for Rv = 3.1 and a gas-to-dust ratio of 200. The column density is simply the integral of density along the 
photon path Nh = J n{r)dl. For a radial path between and r, iVn = K{R~^ — r~^). For a line-of-sight 
path between 6 and Oout = sin~^ [p/Rout), and with impact parameter p, A^h = K{6 — Oout)l P- K is the 
constant of proportionality in the density profile, discussed in §3.3.2, and is in the range 10^^ — 10^'' 
cm~^. 

The extinction in the radial trajectory from the star to the echo is very sensitive to the radius at which 
dust condenses within the CSE, since TVh oc R~^ (provided r ^ Rm)- Using equation (19), and assuming 
p = 1.5 and that dust condenses around = lOOOK, Tc = lOi?* for a 3000K AGB star. This is consistent 
with Habing (1996), who notes that for outflows from cool stars, dust condenses at a radius of 8i?*. At 
this radius and with M = 10~^ Mq yr"-'^, the maximal extinction is Ay = 0.004. For a smaller and hotter 
Cepheid such as RS Pup (i?* ~ 185i?0, ~ 5600K; Deasy 1988), dust wih condense around 60i?*. With 
M = 3 X 10-'^ M0 yr-\ Ay < 0.4, while for M = 10"^ M© yr-\ Ay < A x 10"^. By and large, the 
single-scattering approximation is valid for CSE echoes since the extinction is < 1. 

3.3.5. Background Flux 

Mass-Losing Giants: The entire column of gas containing the echo (delimited with vertical lines) in 
Figure 8 will scatter starlight into the line of sight, contributing the dominant source of background light 
when the echo is imaged at time t. Calculation of this background requires an additional integral of equation 
(7) over the line-of-sight trajectory, also corrected for internal extinction. Consider panel (a) at a later time 
t', say when the variable star is passing through minimum and the echo has left the shell. The echo in the 
column shown will only be detectable in a difference image if the total flux at t (echo and background) can 
be distinguished from the total flux at t' (background only). In the thin-echo approximation, there will never 
be a time when the CSE does not have an echo in it. Instead, one must ask whether the entire column of 
echoes at t is detectable against a background of echoes at a later time t', such as half a period later. To 
implement this, I will use the column of flux at t' as the background flux in equation (16). 

GVs and SNe: I assume that any CSEs these stars might have are optically thin compared to the 
mass-losing giants, and that extinction effects can be ignored. In this case, sources of background flux will 
be the star's quescent light scattered off the surrounding dust, and any unresolved background surface- 
brightness component. The former is estimated as the echo flux diminished by AV for the object (as colors 
change during outburst, this is only an approximation). For the latter, I adopt a background component of 
fly = 18.5 mag arcscc"'^, which is typical of extragalactic disks, in particular the disk of M81 at the location 
of SN 1993 J (Sugerman & Crotts 2002). To scale this to other wavelengths, I use the template spectrum of 
an Sb galaxy from the Kinney-Calzetti catalog within synphot, normalized so the flux in V is unity. 
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3.4. Observability of Light Echoes — Predictions 

Armed with scattering efficiencies S (§3.1), echo constraints (§3.3), and Ught-pulse spectra (§3.2), calcu- 
lating the expected count rate Cobj of light echoes is a simple integration of equation (17) over wavelength for 
a given set of filters and detectors. The resulting exposure times are plotted in Figures 12-16 for a supernova, 
classical nova, TOAD nova, Mira, and Cepheid (respectively). From §3.3.1, "reasonable" exposure times are 
of order Ih for HST and 6h for ground-based imaging, corresponding to logt of 3.5 and 4.3, respectively. 
These have been marked in each figure. 

3.4- 1- Supemovae 

Optical echoes from supernovae are known to be observable, and Figure 12 confirms this for a wide 
variety of echo geometries. Interstellar UV echoes ought to be very difficult to observe, as most lie above 
the exposure-time limits. This is unfortunate since UV echoes are the most effective at constraining the 
grain-size distribution. Since the input spectrum was for SN 1987A 1 day after shock breakout, the UV 
pulse is likely stronger than the adopted spectrum at short wavelengths. This would decrease the exposure 
times in the UV, perhaps making the near-UV echoes observable with long exposures. 

34.2. Novae 

Light echoes from three novae have already been observed (§1), however targeted searches have had 
difficulty making positive detections (van den Bergh 1977; Schaefer 1988). We therefore expect for these 
echoes to be fainter and require longer exposure times. Figure 13a shows that circumstellar echoes may 
be observable across most optical wavebands, but even the closer of the two interstellar echoes [panels (b- 
c)] requires long integration times. If interstellar echoes are to be observed, one is restricted to high dust 
densities, energetic outbursts, and small depths (z < 100 ly). 

3.4.3. Dwarf Novae 

Exposure times for TOAD novae are shown in Figure 14. Circumstellar echoes from accretion-driven 
winds [panel (a)] may be observable in optical wavebands for the brightest TOAD novae, depending on the 
density and thickness of the gas. If such echoes are detected, either in ground-based imaging or with HST, 
we can place much better constraints on the outflow properties. Panel (b) is calculated for the snow-plow of 
the wind into the ISM, and unless the gas densities are unusually high, such echoes, as well as interstellar 
echoes [panel (c)] will most likely not be detected. Given that TOADs are roughly 6 mag brighter than 
ordinary DNe, it is highly unlikely that DNe will produce detectable echoes as well. 

3.4.4. Miras 

Figure 15 shows the exposure times for imaging an echo within the CSE of a Mira. These cool supergiants 
have extremely red spectra and little UV flux, thus we expect their exposure times for the shortest wavebands 
to be quite long. As this is indeed the case, filters shorter than Johnson U are not plotted in this figure. The 
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results are encouraging in that a variety of echoes should be observable. It is of interest however that while 
the template spectrum (Figure llf) has little flux blueward of 6OOOA, echoes from the B and V wavelengths 
may still be observable. This is fortunate since, as shown in §3.1, larger wavelength coverage of an echo 
yields greater constraints on the scattering dust. As a matter of practical importance, these echoes will be 
close (in angular separation) from the central source, which itself is highly luminous and will thus create a 
confusion region. If one can obscure the bright central source (i.e. with a coronagraph) , light echoes can be 
observed within the CSE, provided the dust grain sizes are sufficiently large. 

3.4-5. Cepheids 

Figure 16 shows the exposure times to image variations in nested, thin echoes within the circumstellar 
envelope of a bright Cepheid. As in the previous figure, filters shorter than Johnson U are not plotted. 
For two epochs separated by half a pulsation period, such variations should just be detectable for the high 
mass-loss case [panel (a)] with large platescale detectors. They are effectively undetectable for the lower 
mass-loss rates. Modulations in the total echo signal at a given sightline through the CSE will vary with 
the same frcqiicncy as the pulsating star, which implies that a difference of two images separated by half a 
period will appear as a series of concentric bright and dark rings. This is seen by Sparks (1997) in optical 
images of RS Pup, suggestive of the nested-echo interpretation. Of course, an additional challenge will be 
disentangling the contribution of each nested echo. 

3.5. Caveats 

These results only delineate whether or not echo flux is distinguishable from typical, smooth background 
noise under limited conditions. A second criterium for observability is for the echo to be unambiguously 

distinguishable from background sources. Although in this work I advocate the use of PSF- matched difference 
imaging, I must caution that this technique is optimal only within inherent observational constraints, which 
I briefly enumerate below. (1) To correctly compute the Fourier transform of a PSF, it must have sufficient 
spatial sampling, which translates in practice to a FWIIM> 2.2 pixels. Most HST and some good-seeing 
groimd-based observations (with large detector platescales) produce stellar profiles narrower than this limit, 
necessitating image degradation and introducing subtraction errors. (2) Depending on the CCD used, a star 
that is bright enough to fill a pixel beyond its linear regime can leave subtraction residuals at the core, while 
saturated stars will create residuals within 1-3 FWHM radii. (3) Since point-sources have much higher flux 
than the background, the increased Poisson noise will leave subtraction residuals in the core of any bright 
(but linear) star. (4) Bright point-sources in 775 T observations suffer from large-radius diffraction patterns, 
particularly at short wavelengths. In practice, these PSF-wings are never well subtracted by PSF-matching 
(due to the much lower signal-to- noise in these regions), and PSF-subtraction (i.e. with Tiny Tim model 
PSFs) is only moderately effective at removing this signal at large radii. (5) Many sources of noise can take 
on stellar or extended profiles during data reduction. Hot pixels and cosmic rays mimic undersampled stellar 
profiles following geometric transformation. (6) Charge-transfer efficiency anomalies create extended tails 
around bright stars. (7) Ghosts from bright stars reflecting off optical surfaces produce a wide variety of 
extended structure (e.g. WFPC2 ISR 95-06), including echo-like arcs. 

Given these considerations, an echo should have a sufficiently large angular separation from its source 
to not lie in the confusion region of its PSF, i.e. within a few FWHM. Stars on STIS images, for example. 
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have FWHM~ O'.'l, while ground-based imaging can achieve 0'/5 seeing under extremely-good conditions. 
Let us consider radial separations larger than 0'.'5 to be safe from source confusion. Most echoes imaged to 
date have been rings or arcs of large angular extent (> 20° in position angle). To reliably distinguish an echo 
from a subtraction residual or anomalous background noise, it should be sufficiently well-resolved to have an 
equally-large angular size. At 0'.'5 separation, an echo of arclength 30° will span roughly 6 pixels. Echoes of 
larger arclength should appear as coherant signal over the background and through any other stellar PSFs. 
The chance of observing of an echo is maximized if we require the angular separation between the echo and 
source to be greater than this lower limit of O'.'S, which in turn provides a maximum distance at which the 
echo source can be located. This parameter has been listed in the final column of Table 5, and the value can 
be easily scaled for any other telescope and detector combinations. 

The above arguments imply that the results presented in Figures 12-16 only serve to suggest which 
echoes could be detected. The calculated exposure times do not account for the various noise sources 
discussed, largely because they are much more difficult to characterize in general terms. Stellar subtraction 

residuals are the dominant source of noise, and their presence is a sensitive function of the wavelength 
observed as well as the density and distribution of stars in the field. The readership should fold in these 
considerations with respect to their target fields. Also, I have attempted to generate generic dust geometries 
for producing echoes, but the actual geometry and density can vary substantially from target to target. 
Finally, extremely short exposure times (igxp ^ 1 sec) suggest an echo should be easily detected with an 
integration of typical duration (30-600s), while very long exposure times imply an echo will only be observable 
under very extenuating circumtsances. 



4. Applications 

The primary motivation for building the light echo model presented in this paper was to quantify the 
observability of echoes around a variety of variable objects. Once echoes are detected however, this same 
model is critical in interpreting the data. Given an object's object spectrum, either template or observed, one 
may construct a series of expected spectral distributions against which the actual echoes may be compared, 
to constrain the dust size distribution, composition and gas density. Similarly, if echoes are not detected 
around a given object, this model provides limiting constraints on the above parameters. 

As an application of the techniques described in this work, I discuss the light echoes recently discovered 
around SN 1993 J (Liu et al. 2002, SC02). Both groups report the detection of an echo south-west from 
the SN in WFPC2 imaging taken 8.2 years after core cohapse (UT 19 April 1993; Benson et al. 1994), at a 
radial distance from the SN of roughly 1'.'9 and 120° in angular extent (denoted SW770 in SC02). From the 
echo equation [eq. (1)], this corresponds to a depth z = 770 ly from the SN, assuming a SN distance of 3.63 
Mpc (Preedman et al. 1994). Using PSF-subtraction and PSF-matched difference imaging, SC02 identify 
two additional echo candidates: one located 1'.'15 north-east of the SN (denoted NE260), and a marginal 
detection of SW770 in images taken only 1.8 years after core collapse. The direct and PSF-matched difference 
images in F555W from SC02 are shown in Figures 17a and 17d, respectively. The advantage of difference 
imaging is made most apparent by a comparison of the noise reported by both groups in the F555W image. 
Liu et al. (2002) report a null detection at the position of NE260 with measured noise of 9.2 DN pix"-'^, 
while SC02 measure a flux of 4 DN pix~^ with a background noise of only 1.5 DN pix~^. The reason the 
noise differs between the two groups is addressed in §3.3.1: difl[erence-imaging removes surface brightness 
fluctuations (from non-variable, unresolved stars, nebulosity, etc.) that contribute to the "background" of 
direct photometric techniques. 
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Having constructed a dust model from the the MRN grain-size distribution and the scattering properties 
of Draine & Lee (1984), Liu et aL (2002) present four model fits to their measured surface brightness of SW770 
in F450W, F555W and F814W. Each model overpredicts the blue photometry at ~ 2a, while only the pure- 
graphite model is consistent in the red. The poor fitting in the blue is most likely a result of their overly-dim 
measurement of SW770 in F450W. If one shifts the surface brightness in this filter to within the errors of 
the F555W data (SC02 find roughly equal F450W and F555W fluxes), these models are reasonable fits. Of 
note, they suggest these echoes are explained by dust with density ~ 7 x 10~^° grains cm~^ embedded 
within gas of density nn ^ 3 cm~^. How does the echo model from this work compare? 

In Figure 18, I present model surface brightnesses for a variety of compositions and grain-size distri- 
butions discussed in this paper. In estimating the exposure times in §3.4, I made the assumption that the 
peak-flux spectrum of an object is a reasonable representation of the spectrum over t2, however for de- 
tailed modeling of echoes, this may introduce systematic errors. Integrated fluxes for SN 1993J exist only in 
ground-based filters (Benson et al. 1994; Richmond et al. 1994), however the echoes are measured in WFPC2 
filters, which differ enough from their Johnson and Cousins counterparts that it may be unwise to directly 
intercompare them. I generate the time-integrated flux F{X) as follows. 

Predictions are based on an average of observed spectra from 17 and 24 days after core collapse (Barbon 
et al. 1995, acquired from the SUSPECT supernova spectrum archive). Although taken roughly 3 days 
before and after maximum light, these particular spectra have wider wavelength coverage than that taken 
at maximum light. This spectrum ends around 9500A but the F814W WFPC2 filter extends to llOOOA. As 
an approximate solution, I fill in this gap with the Type la maximum-light template spectrum, renormalized 
so the transition around 9500A is smooth. I caution the reader that this will introduce additional error 
into F814W predictions. This composite spectrum is integrated over the B, V, and Ic response curves, and 
t2 is chosen to match the total ground-based fluence. With t2 — 40 days, this spectrum reproduces the V 
fluence to better than 1%; note that 40 d is also the effective width (/ F{t)dt/ F{tmax)) of the SN light curve 
reported in SC02. The fluences in B and Ic are correct using this value of t2 if the B flux is corrected by 
0.94, and the Ic by 0.77, when integrating over wavelength in equation (17). These corrections are applied 
when computing the predicted WFPC2 fluxes. 

Predicted surface brightnesses are compared to the observed values reported by SC02 for two positions 
along SW770: the central region around PA 220, and the extremal PA 270, at which the echo becomes 
unusually blue. In each panel, the gas density nn has been individually calculated to produce the observed 
F555W measurement for each echo location. I consider pure silicate dust, pure carbonaceous dust, and a 
Galactic composition (equal mixture of both as given by Weingartner & Draine 2001), for the type S, M and 
G grain-size distributions with Ry = 3.1 and 5.5. Toward SN 1993J, Richmond et al. (1994) and Prabhu et 
al. (1995) flnd E{B -V)= 0.08 - 0.32 mag, while values for Ay range from 0.25-1.0 mag (Richmond et al. 
1994; Lewis ct al. 1994), yielding Ry as high as 7.3. For completeness, I also present the model results for 
the MRN size distribution. 

As the model predictions in F814W are the least accurate, these data should carry the least weight 
when considering the goodness of fit. A comparison of the fit to F450W and F555W only shows that all 
three MRN models, and the Type G Ry = 3.1 models for pure silicate and Galactic dust, are the best fits. 
It is of little surprise that the MRN fits closely mimic the Type G Ry = 3.1 fits, since both density models 
were tuned to fit Galactic dust. Including the F814W data as a constraint suggests the need for a highly 
carbonaceous dust distribution. 

We can also judge the goodness of fit by the implied extinction of the predicted dust distributions. With 
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T = 0.11 y, p = 108 ly, and Ap = 11.5 ly, equation (11) suggests the echoing material has a Une of sight 
depth of Az = 47 ly. Using the extinction data in Figure 7, Type G silicate and carbonaceous dust with 
riH = 10 cm~^ will produce Ay=0.40 and 0.33 (respectively) for Ry = 3.1, and 0.49 and 0.20 for Ry = 5.5. 
Thus, Galactic dust will produce respective extinctions of Ay=0.73 and 0.69 for Ry = 3.1 and 5.5. Since 
extinction scales linearly with hydrogen density, purely carbonaceous type G dust would cause roughly 2 
mags of visual extinction. As also noted by Liu et al. (2002), this is an unlikely fit as such optically thick dust 
would obscure background star and galactic light as well. Similarly, the type M and S dust models would 
require unreasonably-high extinctions, and can be safely discarded. With nu = 10 cm~^, and including 
roughly 0.25 mag of Galactic extinction, the Type G Galactic dust models predict Ay ~ 1, consistent with 
the extinction found near the SN by Richmond ct al. (1994). 

It is of little surprise that the above results are inconclusive, since there are many more degrees of 
freedom than constraints. This stresses the need for broad waveband coverage in the interpretation of echo 
signal. In particular, follow-up observations of SN 1993 J should integrate longer and include R and f7-band 
imaging. Despite the shortcomings of this dataset, it illustrates a number of important points. 

Foremost, the observed surface brightnesses arc well reproduced by the light echo model. The dust in 
the disk of M81 (in the vicinity of SN 1993J) appears to follow a grain-size distribution and have a chemical 
composition similar to that of Galactic dust. Gas densities of order 10 cm~^ are required to produce the 
observed echoes, which for Ry = 3.1 translates into dust densities of roughly 10~^ cm~^, in agreement with 
that found by Liu et al. (2002). The color shift in PA between the two echo locations discussed is easily 
explained by a small change in the density of the echoing material. This may be a simpler hypothesis than 
speculating the presence of large compositional gradients of (mostly) pure silicate to pure carbonaceous dust. 
Furthermore, we may also propose an explanation for the absence of the detected echo at other PAs. The 
noise reported by SC02 for their F555W difference image is 1 DN pix^-'^ for a scaled exposure timc^ of 953 
s. This corresponds to a la detection limit of 25 mag arcsec"^, or a hydrogen density of about 2 cm^"^. By 
no means conclusive, this is suggestive that a simple density gradient in the cloud or sheet containing the 
echoing dust can adequately explain the observations. This substructure may be better understood through 
repeated observations and with a more extensive wavelength coverage. 



5. Conclusions 

In this paper, I have presented a model for calculating the brightness of a scattered light echo. I 
find that echoes from supernovae should be the most easily observed, while cataclysmic variables may only 
produce detectable echoes in the innermost circumstellar regions. Post-AGB and Cepheid variables should 
also produce observable echo signal, although the interpretation of such data will be technically challenging. 

These results help explain why previous targeted searches have had such difficulty in directly detecting 
light echoes. For novae, it is fairly likely that echoes were too faint to be observed. For supernovae, the 
case is not as clear. In itself, a genuine null result is interesting, since it places reasonably strong limits 
on the dust density distribution along the line of sight. More likely is that many echoes were lost within 
the surface brightness fluctuations of background unresolved starlight (for SNe in external galaxies) or even 



^Difference images are empirically scaled to stellar flux of the reference image prior to subtraction. In this case, the input 
was a 2000 sec F555W WF4 integration taJcen 2001 June 4, and the reference was a 900 sec F555W PC integration taken 1995 
January 31. A scaling of 953 sec includes the effect of resampling the PC image to WF4 resolution. 
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moderately-crowded stellar fields, typical of observations close to the Galactic plane. 

Consider the effect of simply increasing the platescalc of the detector in otherwise ideal observing 
conditions. As a single pixel images a larger portion of the sky, a light echo will fall on fewer pixels, 
becomming increasingly confused with the background. Although the echo flux per pixel increases (aiding 
detection), so does the background, and point sources become increasingly crowded. The echo SW770 around 
SN 1993J (§4, Fig. 17a), happens to be clearly visible above the background in direct F555W imaging 
(m ~ O'.'l pix"-'^), but is completely lost when the resolution is degraded by a factor of four, as shown in 
Figure 17c. A platescale of 0'.'4 pix~^ is not uncommon by ground-based standards, and for such an echo, 
direct photometric techniques would likely have yielded a null detection. Unlike SW770, echo NE260 is 
only visible in the difference image, since in the direct data, its flux is indistinguishable from background 
fluctuations. Even with the high resolution of HST, direct imaging will miss echo signal. 

Now consider the same exercise using PSF-matchcd difference imaging. Since all sources of constant 
brightness will be removed, the echo still stands out notably from the background when degraded images are 
differenced (Fig. 17e-f). Still at issue is resolving the echo from its source, which we see can not be done for 
NE260 at 0'.'4 pix"-'^. A successful scries of observations must therefore have sufficient resolution to separate 
the echo from its (variable) source, although this does not forcedly require the use of HST. In the case of 
SN 1987A, ground-based imaging was more than sufficient to image interstellar echoes (Xu et al. 1995). For 
ideal difference-image quality, effort should be taken to ensure a stable and well-focused PSF, with individual 
exposures timed to avoid bright star saturation, and total exposure time sufficient to minimize background 
noise. 

It is noteworthy that well-resolved echoes in ground-based observing can be just as easily imaged as from 
space. This is mainly due to the platescale factor of m\ in equation (17). A ground-based integration will 
have > 9 times the flux per pixel compared to the ACS WFC camera, which more than compensates for the 
increased sky background over vacuum observations. Additionally, ground-based observations generally have 
a wider field-of-view, thereby maximizing the volume probed for light echoes. However, increased platescale 
corresponds to diminsished angular resolution, requiring the echo source to be either closer to earth or more 
distant on the sky. Also, seeing greatly diminishes spatial resolution, thereby washing out any substructure 
within the echo. The observer must decide which consideration is more important. 

The question of when to observe an echo remains. Unfortunately, there is no concrete prescription. 
The best image to difference against is one in which there are no echoes, which implies either observing the 
object prior to a burst, or median-combining a long baseline of imaging to attempt to average out any echo 
signal. Immediately following e.g. a nova outburst, the central source is too bright for deep imaging, yet it is 
during this time that echoes sweep through the innermost circumstellar material. Coronagraphy is likely the 
best options at these early times. Long after an outburst, the most easily observed echoes will come from 
forward-scattering geometries, which for large t imply large depth z and angular separation p. However, 
since echo brightness decreases with increasing dust distance and scattering angle, and suffers increasing 
extinction with total optical path length, echoes will become fainter as time progresses. While only a small 
effect for recently discovered SNe, this must be considered for Galactic SNe such as SN 1006. Echoes from 
historic outbursts must be evaluated on a per-case basis with the model presented here. Outside these two 
extremes, imaging should be taken in regular intervals of roughly twice the outburst duration, to minimize 
the overlap of a single echo off the same patch dust in both epochs. With such considerations, light echoes 
should be detected around a much broader list of objects than currently exists. To this end, an observational 
effort to put these ideas into practice is currently underway for the four classes of objects discussed in this 
work. 
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Fig. 1. — Echo geometry used in this paper. Note that the distances are not to scale, and D is much larger 
than z. Adapted from Xu et al. (1994). 
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Fig. 2. — Schematic showing the contribution of dust thickness and finite light-pulse duration on the observed 
width of a light echo. Adapted from Xu et al. (1995). 
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Fig. 3. — Geometric factors in Bsc [cq. (7)]. (a) The geometric function G(p, t, r, Az) [cq. (12) 
(t, r, Az) as follows: solid line (1 y, 3 d, 0.1 ly); dashed (1 y, 0.25 y, 0.1 ly); dot-dashed (10 y, 
dotted (10 y, 0.25 y, 5 ly); dot-dot-dot-dashed (1000 y, 0.25 y, 25 ly). (b) e{p,t) for t=l y (solid 
(dashed line), 1000 y (dotted line). 
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Fig. 4. — Light curves for variable objects. Solid curve - mean light curve below maximum light in B for 
Type ITL SNe. Points — V light curve (shifted by -5 mags) for the Mira x Cyg (from AAVSO International 
Database). The maxima can be well represented by semi-logarithmic parabolae (§3.3.1). 
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Fig. 5. — (a) Scattering efficiency Qsc(A, a) and (b) forward scattering g(A, a) for carbonaceous (heavy lines) 
and silicate (thin lines) grains. Solid lines show values for lOOOA, dashed lines for 2500A, dot-dashed lines 
for 5000A, and dotted lines for 8900A. For clarity, values for the smallest grain sizes are shown at different 
scales, as noted on the ordinates and abscissae. 
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Fig. 6. — Dust grain density distributions for carbonaceous (solid lines) and silicate (dotted lines) dust. 
Thick lines are for the Ry = 3.1, 6c = 6 x 10~^ model, and thin lines are for the Ry = 5.5, 6c = 3 x 10~^ 
model of WDOl. For clarity, and to allow the reader to estimate the mass present in each size range, data 
are plotted in units of a'^ dugr/da. The dot-dashed line shows the corresponding density distribution from 
MRN. 
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Fig. 7. Extinction curves as functions of Ry. The lower curves show the extinction per unit hydrogen 
column density for carbonaceous (solid line) and silicate (dotted line) dust as integrated from equation (18) 
(§2.6). The thick upper line is the total extinction, i.e. the sum of these two components. The thin upper 
line is the extinction from Cardelli et al. (1989), normalized so that Ay = 5.02 x 10~^^ for Ry = 3.1, and 
Ay = 4.77 x 10-22 for Ry = 5.5. 
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Fig. 8. — Schematic diagrams of light echoes within CSEs. (a) A thick-echo regime where only one echo 
passes through the CSE at a time (typical of Miras). (b) A thin-echo regime where many echoes fill the CSE 
at once (typical of Cepheids). See text (§§2.7, 3.3.4 and 3.3.5). 
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Fig. 9. — Integrated scattering function S{\,^) [eq. (8)]. Line thickness designates the density distribution 
model: thick lines are for Rv = 3.1, and thin lines for Rv ~ 5.5 (§2.4). Line style designates dust composi- 
tion: solid lines show carbonaceous dust, and dotted lines show silicates. The first row has been calculated 
for type S dust (§3.1), the second row for type M and the third row for type G. The first column is calculated 
for e = 0°, the second for 6 = 30°, and the third for 6 = 180°. 
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Fig. 10. — Duration of brightening r versus change in brightness AV for the objects listed in Table 4 
Supernovae are marked with stars, classical and recurrent novae with triangles, cataclysmic variables (dwarf 
and TOAD novae) with squares, and mass-losing giant stars with circles. 
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Fig. 11. — Spectra used in this work, with fluxes normalized to a distance of 1 pc. Phoenix synthetic spectra 
for (a) SN 1987A at day 1, and (b) at day 58, (c) a classical nova during the fireball stage, and (d) a classical 
nova during the constant luminosity phase; (e) the M7III star SW Vir (typical of Miras) from the BPGS 
catalog within synphot; (f) the FOIV star ^ Ser from the BPGS catalog, normalized to the F6Ib spectrum of 
HD 8992 (typical of a bright Cepheid). 
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Fig. 12. — Exposure times necessary to image an echo with S/N — 3 (per pixel) for a supernova. Filters 
are plotted at their central wavelengths (Table 2), with STIS filters denoted by crosses, Johnson/Cousins by 
open squares, and ACS filters by filled symbols, with stars for the SBC, squares for the HRC, and triangles 
for the WFC. Geometric parameters are listed at the top of each panel and correspond to those in Table 
5. The dust is assumed to follow a type-G (Galactic) size distribution. Dashed and dotted lines denote 
"reasonable" exposure-time upper limits for HST and ground-based observations, respectively. UV filter 
((A) < 2600 A) exposure times are for the UV burst spectrum, all others are for the optical peak spectrum. 
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Fig. 13. — Same as Figure 12 but for a classical nova. As in Figure 12, UV filter ((A) < 2600A) exposure 
times are for the UV burst spectrum, all others are for the optical peak spectrum. 
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Fig. 14. — Same as Figure 12 but for a TOAD nova. There is no distinction between UV and optical exposure 
times. 
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Fig. 15. — Same as Figure 12 but for a Mira, assuming the thick-echo approximation (§2.7). UV exposure 
times are all longer than the ordinate limits. 
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Fig. 16. — Same as Figure 12 but for a Cepheid, assuming the thin-echo approximation (§2.7). UV exposure 
times arc all longer than the ordinate limits, (a) is for a mass-loss rate of 10~^ Mq yr~^, and (b) is for 
10-9 M0 yr-^ 
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Fig. 17.— F555W images of echoes around SN 1993J from SC02. The SN (at the origin) has been PSF 
subtracted but residuals remain. Panel (a) shows the direct image on WF4 taken in 2001, with resolution 
degraded by two in (b) and again in (c). Panel (d) shows the PSF-matched difference between 2001 and 
1995, and difference images between the input images degraded by two in (e) and again by two in (f). 
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Fig. 18. — Model predictions for the brightness of echoes around SN 1993J for compositions as noted at the 
top of each column. Solid points are surface brigtnesscs as measured by SC02 for the SW770 echo, plotted 
about the central wavelengths of the corresponding WFPC2 filters (F450W, F555W, and F814W). Open 
points are the model predictions. Within each pair, the left point (circles) is for the echo near PA 220, and 
the right point (triangles) is for PA 270. Each model fit is scaled by nu calculated to match the F555W 
datum. These values of nu, as well as the residuals and the grain-size distribution, are noted in each 
panel. 
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Table 1. Mean Widths of Supernova Light Curves 



SN Type 






T (days) 




% of total flux included 




1 


2 


3 4 


5 


1,2 3 4 5 


la 


24 


22 


31 38 


46 


65% 74% 80% 90% 


II-P 


46 


42 


64 120 


160 


76% 65% 95% 95% 


II-L 


26 


28 


39 60 


72 


94% 95% 95% 96% 



^Effective width that Sparks (1994) measured starting at maxi- 
mum hght. 

^Effective width that I measure starting at maximum hght. 
^Effective width that I measure using the entire light curve. 
*Time for Ught curve to rise and fall by 2 mags. 
^2t2 as measured starting from maximum light. 
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Table 2. Characteristics of Filters Used in This Work 



Sysytem 


Filter 




A \ b 

AA 


'T^ C 

peak 


d 

^sky 


STIS 


FUV MAMA 


1381, 


.9 


324.0 


0.045 


3.5 X 10~^ 


STIS 


NUV MAMA 


2309, 


.5 


1237. 


0.031 


0.025 


STIS 


F28X50LP 


7318. 


.7 


2684.9 


0.124 


0.01 


ACS WFC 


F435W 


4321, 


.9 


698.7 


0.372 


0.028 


ACS WFC 


F555W 


5375, 


.9 


847.5 


0.375 


0.055 


ACS WFC 


F606W 


5953, 


.6 


1575. 


0.468 


0.136 


ACS WFC 


F814W 


8092, 


.2 


1541. 


0.439 


0.104 


ACS HRC 


F220W 


2265, 


.6 


442.1 


0.051 


10-4 


ACS HRC 


F330W 


3366, 


,2 


410.7 


0.105 


6 X 10-"^ 


ACS HRC 


F435W 


4321, 


.2 


728.5 


0.223 


4.3 X 10-3 


ACS HRC 


F555W 


5369 


.6 


841.4 


0.239 


9.9 X 10-3 


ACS HRC 


F606W 


5920, 


,1 


1557. 


0.276 


0.024 


ACS HRC 


F814W 


8152, 


.3 


1651. 


0.223 


0.016 


ACS SBC 


F115LP 


1415, 


.0 


354.5 


0.0580 


0.024 


ACS SBC 


F125LP 


1445, 


.3 


334.1 


0.0529 


0.0027 


Johnson 


U 


3673, 


.0 


478.0 


0.050 


0.17 


Johnson 


B 


4519, 


.2 


894.3 


0.058 


0.13 


Johnson 


V 


5556, 


.6 


836.3 


0.115 


0.26 


Johnson 


R 


6784, 


.6 


1559. 


0.137 


0.98 


Cousins 


Ic 


7818, 


.9 


873.8 


0.094 


0.90 


Johnson 


I 


8285, 


.7 


1464. 


0.076 


0.86 



Average wavelength (A) = / T(A)AdA/ / T{X)dX in A. 
^FWHM of filter in A. 
'^Peak throughput of filter and telescope. 

''Sky background in e- s-^ pix-^, for "typical" conditions defined in 
§2.5. 
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Table 3. Characteristics of Detectors Used in This Work 



Instrument 


Platcscalc 


Read Noise 


Gain 


Dark Current 




arcscc pix^^ 


e^ 


e" DN-^ 


e^ s^^ pix^^ 


STIS FUV MAMA 


0.0246 


0. 


1. 


3 X 10-5 


STIS NUV MAMA 


0.0246 


0. 


1. 


2 X 10"^ 


STIS CCD 


0.0507 


4.2 


1. 


2.5 X 10"^ 


ACS WFC 


0.05 


5.0 


1.0 


2 X 10-3 


ACS HRC 


0.027 


4.7 


2.2 


2.5 X 10-3 


ACS SBC 


0.032 


0.0 


1.0 


1.2 X 10-5 


SITE f/lS"* 


0.15 


5.4 


3.3 


4. X 10-3 



Ground-based observations on a generic f/15 2.4m ground-based telescope. See 
i2.5 
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Table 4. Variable Objects with > 2 mag 



Class Example AV (mag) W (days)* Reference 



Supernovae 



SN la 




18 


40 


1 


SN II-P 




17 


140 


1 


SN II-L 




17 


70 


1 


CVs 


CI Nova 


LMC 1991 12 


20-25 


2 


U Sco 


U Sco 


9.2 


16 


3 


RS Oph 


RS Oph 


6.9 


6-8 


3 


U Gem 


U Gem 


6.1 


17.8 


3,4 


Z Cam 


WW Get 


6.4 


9 


3,4 


SS Cyg 


SS Cyg 


3.2 


8 


3,4 


SU Uma 


SW Uma 


8. 


23.1 


3,5 


WZ Sge 


VY Aqr 


9.5 


20 


3,5 


Symbiotic 


Z And 


4.4 


200 


6 


Pulsating Giants 


Mira/LPV 


X Cyg 


8. 


100 


6 


RV Tauri 


R Set 


3.2 


100 


6 


Post-AGB 


CW Leo 


3.8 


220 


6 


Cepheid 


AS Mus 


2.2 


8 


6 


Eruptive Stars 


S Doradus 


r] Car 


9 


2020 


6 



''Duration of outburst or r. 

References. (1) composite light curves from 

Doggett & Branch (1985); (2) Schwarz et al. (2001); 

(3) Catalog and Atlas of Cataclysmic Variables 
(http: / /icanis. stsci.edu/~downes/cvcat/index. html); 

(4) Ak, Ozkan, & Mattci (2002); (5) HowcU et 
al. (1995); (6) General Catalog of Variable Stars 
(http: / /www.sai. msu.su/groups /cluster / gcvs / gcvs/) 
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Table 5. Adopted Geometric Parameters 



Echo"* 


P 


t 




z 


9 


Ap 


^max 




(ly) 


(y) 


(y) (ly) 


(ly) 


n 


(ly) 


(kpc) 




Supernovae 


C 


0.2 


2. 


0.003/0.383 0.05 


-0.99 


168 


2. 


25. 


C 


3.0 


5. 


0.003/0.383 0.10 


-1.6 


118 


0.47 


380 


B 


6.0 


3. 


0.003/0.383 1.0 


4.5 


53. 


0.69 


760 


B 


16. 


5. 


0.003/0.383 2.0 


23.1 


35. 


0.92 


2.0 X 10^ 


I 


60. 


5. 


0.003/0.383 50. 


358 


10. 


4.8 


7.6 X 10^ 


I 


200. 


6. 


0.003/0.383 200. 


3330 


3.42 


8.8 


25. X 10^ 


Novae 


C 


2.5 


0.6 


0.003/0.066 0.5 


4.9 


27. 


0.19 


320 


I 


4.4 


0.19 


0.003/0.066 5.0 


50. 


5. 


0.80 


560 


I 


21.9 


0.95 


0.003/0.066 25.0 


250. 


5. 


1.3 


2.7 X 10"^ 


TOAD Novae 


C 


0.022 


0.0034 


0.05 0.0027 


0.069 


30. 


0.17 


0.3 


B 


0.88 


0.08 


0.05 1. 


5. 


10. 


0.29 


110 


I 


4.3 


0.38 


0.05 2.4 


24. 


10. 


0.30 


10^ 


Miras 


C 


0.0159 


0.141 


0.137 0.097 


0.027 


30. 


0.002 


2.0 


c 


0.0529 


0.151 


0.137 0.136 


0.092 


30. 


0.002 


6.7 


Cepheids 


C 


0.264 


0.0977 


0.0274 N/A& 


N/A 


N/A 


0.05 


33. 



'^C: Circumstellar Echo; B: Contact Discontinuity; I: Interstellar Echo 

'^If two values listed, the first is for the UV burst, the second for the optical maxi- 
mum. For Giant Stars, this column lists Ap in ly. 

■^Maximum distance of source assuming a mimimum echo/source separation of 0'.'5 
using the HST. Sec text (§3.4). 

''For the thin-echo case, the sum of all echoes within the CSE is measured, thus 
Az, z, and 6 are not applicable. 
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